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Abstract 



The issue of the parameterization of small-scale dynamics is addressed 
in the context of passive-scalar turbulence. The basic idea of our strategy 
is to identify dynamical equations for the coarse-grained scalar dynamics 
' starting from closed equations for two-point statistical indicators. With 

the aim of performing a fully-analytical study, the Kraichnan advection 
model is considered. The white-in-time character of the latter model in- 
deed leads to closed equations for the equal-time scalar correlation func- 
tions. The classical closure problem however still arises if a standard 
filtering procedure is applied to those equations in the spirit of the large- 
eddy-simulation strategy. We show both how to perform exact closures 
and how to identify the corresponding coarse-grained scalar evolution. 

j— I . 

^ ! 1 Introduction 

, A common feature characterizing turbulent systems is the presence of space 

■ and time fluctuations 1 . Depending on the strength of some control param- 
eters (e.g. the Reynolds number in hydrodynamic turbulence or the Peclet 
number in scalar turbulence) the number of active degrees of freedom can reach 
values arbitrarily large and eventually diverge in the limit of the so-called fully 
developed turbulence regime pQ. 

.£h [ If, on the one hand, this huge number of active scales characterizing the latter 

■ regime is the ideal framework to investigate 'classical' problem in the realm of 
basic turbulent research (e.g. in relation to global scale invariance and its viola- 
tion in the form of intermittency and anomalous scaling) , on the other hand the 
proliferation of degrees of freedom leaves severe limitations to the deterministic 
description of turbulent fields. To have a quantitative idea, in the turbulent 
atmospheric boundary layer, the degrees of freedom to describe all active scales 
can reach values of the order of, or larger than, 10 27 (10 9 in each spatial direction 
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2 ). This number is obviously far from being handled even by state-of-the-art 
super computers. 

As a matter of fact, in many situations of practical interest (e.g. the description 
of the evolution of a pollutant emitted by sources in the atmospheric boundary 
layer) one is not interested in describing the details (i.e. small-scale dynamics) 
of turbulent fields but, rather, in focusing on their large-scale behaviour. This 
is actually only an apparent simplification. Indeed, the nonlinear character of 
many turbulent systems (e.g. hydrodynamic turbulence) does not permit to 
restrict the attention on the sole (large) scales of interests: small scales enter 
large-scale dynamics in a way that the latter is not closed. This is the well-known 
closure problem j^j. It is worth stressing that closure problems also appear in 
the absence of nonlinearity. The best example is provided by passive-scalar tur- 
bulence, where the closure problem arises from the advection term (i.e. a term 
involving a multiplicative structure between velocity and scalar). 
The present situation, in relation to the long-standing closure problem, is that 
a unifying theory predicting how to close the basic equations in the large scales 
is not available. Readers can refer to 3 for a modern review on the closure 
problem oriented to the so-called large-eddy simulation (LES) strategy. 
The situation is even more difficult when different turbulent systems are sup- 
posed to interact with each other. One of the best-known examples is the 
temperature field in a turbulent flow, which, depending on the role of buoyancy, 
may behave as an active or a passive scalar. If very few exact results on the 
closure problem are known for non interacting turbulent systems, the situation 
is even worst in the presence of interaction. 

The aim of the present paper is to provide some exact results for a particular 
class of turbulent systems: the passive-scalar turbulence. This problem is per se 
interesting in connection, e.g., to numerical studies of pollutant concentration 
in the atmosphere, but it might also give useful suggestions on how to generalize 
the results to active scalar fields and, hopefully, to hydrodynamic turbulence. 
This latter point is however, at the present stage, far from being achieved. 

The paper is organized as follow. In S|2]we recall the basic features of passive- 
scalar turbulence. In we introduce the large-eddy simulation technique by 
defining the filtering process, the coarse-grained fields and their subgrid coun- 
terparts. The general unclosed character of the equations does not prevent the 
existence of situations in which exact closures are possible, an example of which 
is shown in [JU Starting from fJSwc focus our attention on the Kraichnan ad- 
vection model and, in ^ we analyse the filtered equation. provides a series 
of possible closures, emerging with different degrees of approximation or in dif- 
ferent contexts. Conclusions and perspectives follow in SjSJ Appendices flAl and 
flBl are devoted to provide some exact analytical expressions and to study the 
limit cases of purely diffusive or smooth flows, respectively. 
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2 Passive scalar turbulence 



The basic equation governing the dynamics of a passive scalar field 8(x, t) is the 
well-known advection-diffusion forced equation 

d t 0(x, t) + v(x, t) ■ d8(x, t) = n Q d 2 9(x, t) + f(x, t) . (1) 

The advecting velocity field v(x, t) is assumed incompressible: d ■ v = 0. Scalar 
fluctuations are injected into the system at a large scale L by the forcing term 
f(x,t), acting as an external source. Scalar dissipation takes place at small 
scales, of order 77, and is accounted for by the molecular diffusivity k . 
A strong analogy with the Navier-Stokes (NS) turbulence holds: the number of 
active spatial degrees of freedom 2 rapidly increases with the Peclet number Pe 
(the analogous of the Reynolds number) and scales as {L/ r q) d ~ Pe 3d ^ 4 , where 
d is the space dimension (> 2). The advantage of dealing with (JIJ instead of 
NS consists both in the linearity of the equation and in the scalar character of 
the unknown field, and finally in the locality of the physical-space description 
(i.e. the analogous of the pressure field is absent). Nevertheless, the importance 
of scalar turbulence is underlined by the considerable progress that has been 
achieved in the last few years in this context 0|. 

3 Large-eddy simulation 

3.1 Definition of the filtering process 

Coarse-grained fields (denoted with a tilde) are obtained from the original, fully- 
resolved fields through a convolution with a low-pass isotropic filter Pi of char- 
acteristic length ?, lying in the inertial range of scales: 77 -C I <C L. In particular, 
we have 

§(x,t) = Jd d yP l (x-y)6(y,t) = Jd d sP l {s)e{x + s,t) = {Pi*6)(x,t) , 

and similarly for v = Pi* v and / = Pi* f. The filtering process thus defines 
a linear operator which commutates with any partial derivative, because of the 
structure of the convolution kernel, but which does not factorize when acting 
on products (the filtered of a product is not the product of the filtered) . 
Small-scale fluctuations, denoted with a star, are defined as: 

0* = 6-6, v*=v-v, /* = /-/. (2) 

3.2 The problem of closure in the large scale 

The large-eddy simulation (LES) strategy is carried out by convolving 0) with 
the filter P/, in order to obtain an equation for the coarse-grained fields: 

d t (i + : v^d9 = K d 2 6 + f . (3) 
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Unfortunately, this is not a closed equation in the large scales, as the multi- 
plicative (in the filtering operation) term still involves a product between fully- 
resolved fields. We thus rearrange it in the form 



d t e + v-d6 = K Q d 2 9 + f-d-r e , 



(4) 



where the subgrid scalar flux Tg is given by: 



Tg = W 



9 



V' 
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(5) 



The aim of LES closures |3J is to express d-Tg in terms of coarse-grained fields, in 
order to get a closed equation describing the large-scale dynamics autonomously. 
Once this purpose is accomplished, Q can be integrated numerically (see, e.g., 
0) on a mesh of spacing / instead of 77, as it would be necessary for the original 
equation This less-demanding integration means a huge gain in memory 
and CPU time requirements and represents the essential advantage of the LES 
strategy. 

From a general point of view, the perfect closure (i.e. having the 'true' Tg) is 
able to correctly represent all the obscrvables built by filtering the 'true' field 6. 
On the contrary, only some observables, and some of them in principle better 
than others, can be correctly described by any empirical closure. It is worth 
noticing that, which observable is properly reproduced, can be assessed only a 
posteriori (e.g. by comparing LES predictions against experiments). 
Unluckily, no general closed expression for d-Tg in terms of 9 and v is available: 
this is a clear indication of the strong coupling between all scales which is typical 
in turbulence. A remarkable exception is provided by the case where there is a 
marked scale separation between velocity and scalar length and time scales. It 
is then possible to show [7| |S] that the effect of unresolved scales is just the 
rcnormalization of the molecular diffusion coefficient k,q to an enhanced eddy 
diffusivity K c g (generally speaking, an eddy-diffusivity tensorial field). General 
expressions for the eddy diffusivity as a function of the flow properties do not 
exist, and in most cases K c ff can be determined only numerically. 
Here, our aim is to consider the challenging situation where there is no scale 
separation 9 between velocity and scalar and to explore, in such a context, the 
existence of effective equations for 9. 

3.3 Structure of small-scale contributions 

An alternative expression of d-Tg can be obtained, for future purpose, plugging 
in © the decomposition J3J). We have: 



d-rg=C + C + TZ 



(6) 



where 




(Leonard-like term) 
(cross term) 
(Reynolds-like term) 



(7) 



4 



4 An example of exact closure 



We now specialize to stochastic velocities and forcings, i.e. to the case in which 
both v and / are fields with assigned statistical properties; we denote with 
brackets the average over their statistical distribution: it is then clear that 
this ensemble average commutates with the spatial average represented by the 
filtering process. 

Also in this framework, a renowned example of exact closure exists. Indeed, for 
times t larger than the largest velocity time scale, the mean field (9) experiences 
the cumulative effect of velocity via an eddy-diffusivity coefficient: 

d t (6) = n tot d 2 (9) + (/) . (8) 

The expression for «tot (here supposed to be isotropic, as actually is in the 
presence of isotropic velocity fields) follows from the well-known Taylor formula: 

1 f°° 

Having accounted for the advective term through a linear diffusive one, in this 
particular case the nonlinearity problem faced while passing from equation Q 
to yj does not exist any more and the same eddy-diffusivity equation must then 
hold also for the averaged coarse-grained field (9) . By virtue of linearity, from 
iJSJ we have: 

d t (9) = K tot d 2 (9) + (/) . (10) 

We now look for an eddy-diffusivity-type closure in the equation for 9 such that, 
starting from it, (|1U|) is recovered. We ask, in other words, that the closure is 
able to reproduce the averaged, long-time behaviour of 9. The searched equation 
is 

d t 9 + i-d9 = K cS d 2 9 + f, (11) 
where n c s has to be determined. We can use Taylor's formula again to obtain 

1 f°° 

Ktot = Kcff + ^y dr (v(t) ■ v(0)) . (12) 
A simple comparison between equations © and (|12|l yields: 

«c ff = «o + Ydj™ dT [{v{t) ■ v{0)) - ( * (r) ' * (0))] • (13) 

5 Kraichnan advection model 
5.1 Properties of velocity and forcing 

We further restrict our attention to the well-known Kraichnan advection model 
12313) m which both velocity and forcing are Gaussian, white-in-time and zero- 
mean random fields, statistically stationary, homogeneous and isotropic (for an 
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extension to inhomogeneous or anisotropic situations see, e.g., |r2IH5] 1. 
The statistics of v is thus fully determined by its second-order increments, 

{[v^r,t)-v^O,0)}[v v (r,t)-v v (O,0)})=26(t)D^J(r) , 

following an inertial-range spatial power law: 



DW (r) = D r* [(d + £ - 1)8^ ^ j • (14) 

The assumption of (5-correlation in time is of course far from reality, but it 
has the remarkable advantage of leading to closed equations for the equal-time 
correlation functions = (0(xi,t) ■ ■ ■ 9(x n , t)) of any order n (see, e.g., (15)1. 
The parameter f, lying in the open interval (0,2), governs the roughness of 
the velocity field, whose Holder exponent is £/2. Due to the lack of temporal 
memory of the flow, the Kolmogorov value ^1] is £ = 4/3; the limit cases £ = 
and £ = 2 will be studied in the appendix. The diffusive scale rj, at which 
diffusive and advective effects are comparable, is defined by the relation: 

A convenient choice for / is to assume a step-function form for its two-point 
correlation (/(r, t)f(0, 0)) = 5(t)F L (r), i.e. F L {r) = F Q(L - r) (the details 
of its behaviour around the large scale L are not relevant, but in this case the 
inertial range has an upper limit r < L instead of r <C L: see, e.g., |15|1. 

5.2 Two-point correlation function of the fully-resolved 
passive scalar 

In this framework the mean value (6) shows trivial dynamics, while the equation 
for the two-point equal-time correlation function C\ can be derived analyti- 
cally from (JTJ exploiting Furutsu-Novikov-Donsker's (FND) rule for Gaussian 
integration by parts [T71 IT31 IT31 IT] : 

d t ci 8) = flWfl^^' + 2n Q d 2 C { i ] + F L . (16) 

The stationary version of 116fl . in which the left-hand side vanishes and is 
only a function of the distance r between the two points (r = x — x'), can be 
simply solved by splitting the integration interval and gives |2(J) : 

for r -C t] 

C^\r) ~ \ c - kr 2 -^ for r? -C r < L (17) 




where 



for r > L 



F L 2 -t , £F oV 2 -Z 



(d - 1)(2 - 0(d + £ - 2)A) ' 2d(d - 1)(2 - ODo 

d{d-l){2-€)D Q ' 2rfi ' d + £-2 
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In the limit of small viscosity, which will be assumed throughout the paper, the 
merged-point value of the correlation is given by (9 2 ) ~ c and the stationary 
second-order structure function turns out to be a pure power law in the inertial 
range: 

sf (r) = ([6(r, t) - 9(0, t)} 2 ) = 2(9 2 ) - 2C^\r) = 2kr 2 ^ . (18) 

The exponent 2 — £ coincides with the predictions based on dimensional argu- 
ments and becomes 2/3 for £ = 4/3, according to the Kolmogorov-Obukhov- 
Corrsin scaling |21|. 

The equation for the steady-state dissipation, which we write for future purpose, 
arises from (|16|l evaluated at merged points: 

2n ((d9) 2 ) =F . (19) 

Exact solutions, i.e. far from the perturbative limits £ — > 22 , £ — > 2 123 and 
d — > oo |24j . are not available for higher-order correlations. 

5.3 Two-point correlation function of the coarse-grained 
passive scalar 

To provide a benchmark for the various closures, we first evaluate the stationary 
coarse-grained correlation function by its definition: 

C<f(r) ee (§(x,t)9(x',t)) 

= Jd d yJd d y'P l (x-y)P l (x'-y r )(e(y,t)9(y / ,t)) (20) 

= Jd d sJd d s'Pi(s)Pi( S ')Ci e \\r + s + S '\) . 

Its exact value is reported in the appendix for d = 3 and a top-hat spherical 
filter, Pi(s) = 30(1 — s)/4tc1 3 (from now on, unless explicitly stated, we will 
confine ourselves to this situation). Here we only need to express its Taylor's 
expansion in the parameter l/r, with r lying in the inertial range: 



1+1(2-0(3-0 (21) 

3 _ „ / a 4 ^ /r 6 ' 



, i7 -^-l)(2- 00-0 +0^ 

Clearly, as the separation r increases and becomes much greater than the filter 
scale I, the unfiltered result is recovered. On the other side, a lower limit for 
the physical consistence of the latter expansion can be intuitively identified in 
r > 21, because for smaller separations the two integration domains in l|20|l 
would partially overlap. Expression (|2"TJ> represents, therefore, the best result 
that can be achieved by means of a closure. 
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For a 3-D Gaussian filter, P/(s) = (2nl 2 ) 



-3/2, 



I 21 , the coefficient of the (I ft 



term in JUJ is (2 - £)(3 - 0- For a 2-D top-hat filter, P/(s) = 6(7 - s)/7t/ 2 , it 
becomes (2 - £) 2 /4. 



5.4 Properties of the filtered velocity and forcing 

At this stage, it is also convenient to analyse the behaviour of the filtered ve- 
locity and forcing. By definition, both v and / are Gaussian, white-in-time and 
zero-mean random fields, statistically stationary, homogeneous and isotropic. 
The large-scale character of the forcing reflects in the fact that its two-point 
correlation (f(r, t)f(0, 0)) = (5(i)^-i(r) keeps the same spatial step-like form 

T L {r) = Jd d s Jd d s' fi{a)Pl(af)F L (]r + s + s'\) = F L (r) for r £ I , 

where the interval I = [L — 21, L + 21] becomes negligible for I <C L. 
Second-order increments of the coarse-grained velocity, 



[» M (r,t) -«„((>, 0)][e„(r,t) - v u (0,0)]) = 26(t)D$(i 



are given by 



d d s / d d s' Pt( 8 )Pi(8') [l>W(r + « + «')" + (22) 



(23) 



The exact expressions of the coefficients A(r) and B(r) are quite cumbersome 
(see appendix); a much more useful and meaningful quantity is provided by the 
following power-series expansion in l/r: 



D$(r) = D rt{ f(2 + ^-?- 



l + O 



(4 + 0(6 + 
2«48£> / ? 



O 



(4 + 0(6 + 



(24) 



The latter will be extensively used in the following. 

Similar results also hold for both the 3-D Gaussian and the 2-D top-hat filter, 
but with different numerical coefficients. 



6 Analysis of the filtered equation 

The first step, in order to find closed equations for the large scales, consists in 
deriving the exact equations for the two-point correlation function of the filtered 
field. With the same procedure used when passing from Q to l(T5|l. starting 
from equation (J2J we can write: 

d t (0(x,t)§(x\t)) +2(9(x,t)(v-60)(x',t)) = 
= 2n d 2 (6(x,t)6(x',t)) + T L (\x - x'\) - 2(§(x,t)d ■ T 8 (x',t)) . (25) 
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This is the starting point for our systematic procedure to construct closure 
approximations. The second term on the left-hand side depends on v and 9 
only, but we cannot transform it (yet) into a contribution with the structure 
of the first term on the right-hand side of ijlUj) because, at this stage, we are 
not able to apply FND's rule on it, as we do not know the functional derivative 
of 9 with respect to v explicitly. The last term on the right-hand side of l(2Sjl. 
which is not expressed as a function of coarse-grained fields, is the 'disturbing' 
quantity: our purpose is therefore to find its perturbative expansion in I Jr. 
It is not difficult (although quite lengthy) to prove that: 

(iMe . T .c, () , . ( iV - (L)' + o . (26) 



(4 + 0(6 + W 30 

More specifically, exploiting decomposition © and definitions (7J , we can show 
that l|26|l consists of: 



(C(x,t)9(x',t)} = O 
(C(x,t)6(x',t)) = 



(3 -Q(O + 10£ + 24-2^24 )f ( I 
3(4 + 0(6 + 



30 \r J \r , 
(3-O(C 2 + 10C + 24-2«12)F fl\t ^fl^ 2+5 



W*>WV,t)) = ± 3( 4 + 0(6 + [r) + ° 

(27) 

The first line of (|27|l tells us that the Leonard-type term does not contribute, 
at the lowest two orders, to the equation for the two-point correlation of the 
coarse-grained scalar. Since our closures are derived from this equation up to 
the second order, it follows that the Leonard-type term will not contribute to 
small-scale parameterizations. This fact is not a consequence of the Kraichnan 
advection model but rather seems to hold for general advection models. For 
standard closure models based on single-point quantities, the contribution from 
the Leonard stress in the parameterizations is, generally speaking, non-zero. 
Moreover, it is easy to see that the sum of the last two equations in l(27|) exactly 
coincides with 1)26(1 . in spite of the further convolution required in © on the 
cross and Reynolds-like terms. This is in accordance with the following result, 
that can be obtained after simple algebra: 

/ V- ... ^ . s , , . /M — 



(C(x, t)9(x', t)) +0 - = (C(x, t)6(x', t)) = (C(x, t)0(x', t)) + O 

\r J \r 

The same result holds replacing C with 1Z and has been obtained exploiting a 
further identity, holding for any couple of fields j\ and f -2 in the presence of 
statistical homogeneity: (/i(r,*)/ 3 (0,0)) = (/i(r,t)/ a (0,0)). 
It is also interesting to study when the cross term and the Reynolds-like one 
dominate over each other. At the lowest order, (l/r)*, we should thus compare 
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the absolute values of the two numerical coefficients (functions of £) appearing 
on the second and on the last line of (|27|l : the result is that the former prevails 
for £ > £o) w hh £ — 0.92, and the latter is dominant when £ does not exceed 
this critical value £o. Furthermore, one can easily prove that the two terms 
composing C Q) give the same contribution to the parameterization, each being 
half of the numerical coefficient appearing on the second line of l|27[l . 
Once again, expressions with the same structure of (|26() and (|27() also hold for 
the 3-D Gaussian filter, with different numerical coefficients. On the contrary, 
for the 2-D top-hat filter, it turns out that the second-order coefficient in <|2lfl) 
vanishes. 



7 Examples of analytical closures 

7.1 Importance of small-scale contributions 

Let us now focus on consequences of relation (|26[) : it is immediate to realize 
that neglecting small-scales effects completely, i.e. assuming d ■ tq = 0, makes 
equation 1|25|) unbalanced at order (l/r)^. In other words, if one assumes a 
closure of the kind 

d t § + i> ■ 86 = K d 2 9 + f , 

starting from it one would obtain a non-analytical expansion for the coarse- 
grained scalar correlation, C^\r) = c — kr 2 ~^[l + 0(7/r)^], clearly in contrast 
with the exact result 121fl . 



7.2 Constant-eddy-diffusivity closure 

The first issue thus consists in finding a way to take order {l/rp into account 
properly. For this purpose it is sufficient to notice that the diffusive contribution 
in l|25|) , or equivalently (at this order of approximation) in (|16|) , turns out to be 
proportional to r - ^, precisely 

It is then clear that an effective-diffusivity term, like the one proposed in 
must be able to correctly reproduce the lowest-order contribution of More 
specifically, writing K c ff = kq + k±, the equation 

d t + v-d6 = n cff d 2 + f (29) 

is balanced at order (l/r)^ if 

2«24 , , , 

- = (4 + 0(6 + ^ (30) 

(to prove this, it is sufficient to replace Kg with n\ in Q28H and to compare the 
result with (|26|l ). It is worth noticing that expression H^OJI also follows from the 
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integral in (|13[1 . which substantially amounts to compute the difference between 
total and large-scale kinetic energies in the presence of 5-correlated flows. This 
is in accordance with equation (|24[) . from which we deduce 



«1<V = lim . (31) 



Two more remarks about this kind of closure, which we will call 'constant eddy 
difTusivity' (CED), emerge from the consideration that the fraction (function 
of £) in l|30|) always stands between 1 and 2. First, reminding that from l|15fl 
we have ko — D^rf , a remarkable increase in the transport coefficient is found 
{I r\ => K e s — Ki Ko): this is a typical effect in turbulence. Second, an 
effective dissipative scale comparable to the filtering length I arises for 9: the 
analogy is evident between the roles played by I for 9 and by the molecular 
dissipative scale r\ for the original field 9. 

The equation for the coarse-grained scalar correlation arising from l|29|) has 
exactly the same structure of the fully-resolved corresponding l(l"6|) : 

dtC? = D^d^C® + 2 Kcfi d 2 cf ) + T L . (32) 
The stationary solution of 132(1 in the inertial range is: 



i n\ 2 . /r ■ 



l + -(2-0(3 + 0(-l +0 



(33) 



A comparison with the exact result (|21|l shows that the CED closure is able to 
capture the correct order of the deviation from the fully-resolved scalar corre- 
lation l|17fl . i.e. (l/r) 2 , but with a wrong coefficient. It can be proved that the 
maximum error takes place at £ = 1. 

From (1M2|) we derive the following equation for the steady-state dissipation: 

2 Kcfi ((d9) 2 ) = F . (34) 

Owing to the fact that K c ff ^S> no, the comparison of l|19l) with l|34|l tells us that 
the average of the square gradient of the scalar is much smaller for the filtered 
field than for the original one, as one would expect. 

For the 3-D Gaussian filter or the 2-D top-hat filter, expression (|30|l becomes 

2 1+ ? /3 + A 

or 



r ( 2 + ^ . Do fi 



2 r(2 + £/2)r(3 + e/2) 

respectively, T being Euler's function. Consequently, the ratio ki/DqI^ stands 
between 1 and 10 in the former case and between 1/2 and 1 in the latter. 
Plugging these values of k,\ in equation l|29|l or (|32|l and computing the coarse- 
grained correlation, for the 3-D Gaussian filter we obtain a result similar to (|33|) , 
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without the 1/5 factor in the (l/r) 2 term. However, for the 2-D top-hat filter, 
the second-order coefficient is (2 — £) 2 /4 and exactly coincides with its 'true' 
value, i.e. we get the remarkable result that the error in the approximation is 
automatically pushed at higher orders. 

7.3 Improved closure 

Our aim is now to improve the constant-eddy-diffusivity closure, exact at order 
(l/rp, by introducing a new closure which is accurate up to order (l/r) 2 |25|. 
At this stage, we have no hint of how to implement this closure, differently from 
what happens with the intuitive emergence of the eddy diffusivity; a trivial Tay- 
lor expansion on the turbulent fields would actually prove wrong. 
However, we know that we have to reproduce a term proportional to I 2 in the 
equation for the coarse-grained scalar correlation: it is then reasonable to add, 
on the right-hand side of the equation for a new contribution proportional to 
some power of I. The minimal guess could be represented by the addition of a 
term linear in I, which would be able to generate a quadratic correction in the 
equation for the correlation when applying FND's rule. This guess is however 
ruled out by symmetry considerations, because we would need to introduce some 
additional field (with the dimensions of the vorticity) which cannot appear in 
the equation for 6 for parity conservation. 

The next possibility is thus to add a term quadratic in I, in which the coarse- 
grained fields must appear in the tensorial form vdO. Dimensional consider- 
ations then require the presence of a square length at denominator, but no 
scales other than the filter width can appear, because we are dealing with a 
single-point equation and neither L nor 77 are relevant. A second derivative is 
thus required. If one completely neglects higher orders in I (it is important to 
underline that this ansatz is not trivial at all, because they would give rise to 
spurious contributions at lower orders), then the searched equation must have 
the following form: 

d t 9 + v ■ 86 = K cS d 2 9 + f + I 2 (ad 2 v ■ d~6 + (3dv : dd§ + jv ■ dd 2 ^ . (35) 

The coefficients a, (3 and 7 can be uniquely found by imposing the correct 
description of order (l/r) and, at the same time, the vanishing of any spurious 
modification to order (l/r)^, which has already been captured by the CED 
closure. In other words, we ask that the value of n c s previously found remains 
unchanged in 1|35|) (alternative conditions will be shown later in this section). 
With these hypotheses, we have a — and (3 = 7 = — £/15, so that (|35|) can be 
written as: 

d t e + v-ae = n cS d 2 e + f- ^-i^d^dj) (36) 

to 

(for a comparison with the corresponding nonlinear closure in NS turbulence 
see, e.g., [3^1 E7]). Note, however, that (|36|l does not coincide with equation 
(3.16) reported in |2H]. The reason of such a discrepancy has been found in a 
wrong identification done in to infer the closure (see, e.g., |28p. 
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It is worth noticing that this closure, which has been obtained through Eulcrian 
considerations, has no Lagrangian counterpart, differently from CED. Indeed, 
an expansion in the spirit of Kramers-Moyal would yield a term with the same 
structure of the one we have introduced on the right-hand side of iJSSJ, but 
it turns out that our triplet of coefficients a, f3, 7 appearing in l|3tj[> does not 
satisfy the constraints imposed by Pawula's theorem |291 13TJ] . This fact is also 
related to the breaking of Galilean invariance in (|36|l . which has been obtained 
in the frame of reference where the velocity is zero-mean and is not exportable. 
Moreover, the universality of K e ff is no more present in the triplet, whose value 
has been determined exploiting the explicit solution for c!f ^ . 
The equation for the two-point equal-time correlation function arising from l|36[l 
has now a different structure: 



+ 



where Vyx^ satisfies (w M (0, t)d 1 d\v v (0, 0)) = 5(t)Vyx^ v . 
The stationary solution of i|37|) in the inertial range is: 



1 „ /A 2 /P 2+ « 



1 + -(2-0(8-0 - +0 - 



(38) 



Comparing (|38|l with (|21() we conclude that this new kind of closure is able 
to reproduce the exact structure of the coarse-grained scalar correlation up to 
the second order with the correct value of the coefficient, differently from what 
happens in (|3^|) . We also notice that the error is now 0(l/r) 2+ ^ instead of 
0(1 /r) 4 : in order to balance also this contribution in the proper way, one would 
have to add other terms on the right-hand side of (|36|l , paying attention to take 
unmodified the lower orders which have already been captured. 
A last comment on ll.itit is worthwhile. It turns out that the same dissipation 
equation (|34ll generated from CED still holds. This is in accordance with the 
flux- like structure of the last term on the right-hand side of <|5fi|l . which gives 
no contribution to the equation for the dissipation. Exploiting homogeneity (H) 
and incompressibility (I), we have indeed: 

= -((dj)d^dj)) I ((^cU^cU) . 

Comparing the third and the last member of this chain of equalities, we deduce 
that all terms must vanish; in particular, also the first one is thus zero, which 
proves our assertion. 

For the 3-D Gaussian filter, the correct values of the triplet are a = and 
P = 7 = — £/3. The closed equation for 6 has thus the same structure of l|3^|l 
and leads to a coarse-grained correlation expressed by (|38(l but without the 1/5 
factor near (l/r) 1 , which is exact up to the second-order. The same accuracy is 
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obtained for the 2-D top- hat filter with a = (3 = 7 = 0, coherently with what 
stated at the end of [Hand of iZH 

Before proceeding, we would like to stress the conditions we specified when 
(uniquely) fixing the values of a, and 7 in equation l|35|) in order to obtain (|3(jl) : 
1) correct description of 0(1 /r) 2 , 2) complete neglection of terms containing 
powers of I higher than 2, and 3) absence of any modification of 0(7 /r) ? i.e. of 

toff- 

7.4 An alternative closure 

We are now going to show an example of what happens if the possibility of 
taking filters of (slightly) different length between scalar and velocity is allowed. 
In particular, if we define a new coarse-grained velocity as v = P a i * v (with 
a = ~ 6/(3 + £), i.e. slightly less than unity), then we assert that an 
equation with the same structure as in i|29|) . 

d t 9 + v-d§ = n cS d 2 9 + f , 

is able to reproduce for the correlation the same degree of approximation of 
(|38l) . with an error pushed again at 0(Z/r) 4 and without any need of introducing 
additive terms like on the right-hand side of 135fl . The same expression (|13l) for 
the eddy diffusivity still holds, with v replaced by v\ this amounts to perform 
the same replacing in expression (|31|l for k\ — K c g — kq, whose numerical value 
is now times the value given by (|30J) . in accordance with the substitution 
/ 1—* al. 

7.5 Absence of eddy diffusivity 

We now present some other possible closures, which emerge if one pursues aims 
different from those specified previously. 

From another point of view, the condition 3) fixed at the end of ^7.31 could 
be modified into: 3bis) modification of 0(1 /rp such that K e s is restored to its 
original value kq. In other words, we would now like to pass from equation Q 
directly to 

d t 9 + v-dd = K Q d 2 9 + f + I 2 (a'd 2 v ■ 80 + {3'dv : 889 + j'v ■ 8d 2 9^j , (39) 

imposing that the contribution in parentheses on the right-hand side of Ij39(l is 
able to correctly balance both 0(l/r)^ and 0(l/r) 2 . This goal can be accom- 
plished for 

, g + 2 -g 4 - 9g 3 - 18g 2 + 120£ + 240 

a ~ £(3 + 6(6 + 6 ' P ~ 15^(3 + 6(6 + 

, _ -£ 4 - 9<f - 18<f + 60£ + 120 
7 " 15^(3 + 6(6 + 6 ' 

It is worth noticing that, in this way, the term proportional to I 2 in (|39|l has 
no more an overall flux- like structure, differently from what happened in 1)36(1 . 
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Consequently, it gives rise to a non-zero contribution in the equation for the 
dissipation. 

Conditions 3) and 3bis) represent two completely different points of view: in 
the first one the diffusive term captures 0(l/rp completely, in the second one it 
gives no contribution. Of course, there exists an infinite range of intermediate 
possibilities if one considers both the renormalized diffusivity and the triplet as 
unknowns. 

7.6 Absence of second-order correction 

Suppose now to be interested to measure 'pure' inertial-range scaling behaviour 
by means of a LES strategy. More specifically, referring to (|17|l . our aim here 
is to reproduce the asymptotic behaviour C^^r) ~ r 2- ^ even in the presence 
of finite-size effects, i.e. induced by the filter cut-off I, which would imply the 
additive corrections in powers of l/r shown in 1|21|) . In plain words, we renounce 
to describe the actual coarse-grained correlation function (|21|) . with the aim of 
isolating the scaling behaviour which would be observed at higher-Reynolds- 
number (and thus higher-resolution) numerical simulations. 
Assuming the conditions 2) and 3) mentioned at the end of H7.3I this amounts 
to modify condition 1) into: Ibis) absence of second-order correction in the 
coarse-grained correlation. This aim is accomplished simply assuming a = 
and /3 = 7 = -(3 + 0/30. 

It can be shown indeed that, starting from the equation 

die + v ■ do = Kctf d 2 e + f- ^Apd^^dJ) , 

the coarse-grained correlation turns out to be: 

Cf\r)=c-kr 2 -t 

The deviation from the fully-resolved correlation is pushed at a higher order: 
we have thus obtained a better recovery of the original field in its two-point 
correlation function. 

8 Conclusions 

The main aim of the present paper was to provide a systematic procedure to 
derive closed equations for a coarse-grained passive scalar field. Such equa- 
tions have been shown in fQand, starting from each of them, we computed the 
corresponding two-point equal-time scalar correlation function to which they 
give rise, in the stationary state of the Kraichnan advection model. Specifi- 
cally, in H7.2I we performed an analytical derivation of the 'standard' constant- 
eddy-diffusivity parameterization, which also arises in several other contexts 
but which is usually obtained empirically. Our procedure allowed us not only 
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to quantify the error intrinsic in the parameterization, by constantly keeping it 
under control, but also to refine this closure by introducing a new term ( jjEH in 
the equation for 6, thus obtaining a dynamical parameterization which can be 
seen as the passive-scalar counterpart of the nonlinear closure in NS turbulence. 
In this way, we were able to correctly describe both the asymptotic behaviour 
and the scales of motion r in the vicinity of the filter cut-off I (a limit r > 21 can 
roughly be identified for the validity of this description), as this closure captures 
the 0(l/r) 2 correction and pushes the error at higher orders. A natural question 
then arises: is it possible to generalize the aforementioned procedure, in order 
to correctly describe also higher orders? This problem is twofold, according to 
the meaning of the word 'order'. 

First, always focusing on the two-point correlation, one would like to be able 
to consider separations r closer and closer to the filtering length I (or, cquiva- 
lently, to use a larger, numerically less-demanding, filtering length, to describe 
the same separation): this could be accomplished by balancing equation l|25JI 
also at some further orders, e.g. 0(//r) 2+ ^ and 0(l/r) 4 . This task looks fea- 
sible by adding another term on the right-hand side of (|36f) . proportional to l A 
and whose tensorial structure could be determined by symmetry considerations. 
However, finding the exact numerical coefficients (the counterpart of the triplet 
a, /3, 7) would imply cumbersome calculations, and — in view of what stated 
before equation (|35|l and in condition 2) at the end of H7.3I — spurious effects 
at lower orders may appear. 

A second, and probably much more meaningful, generalization of the present 
work would be represented by the extension to higher-order (always equal-time) 
correlation functions, e.g. C4. Our analysis has been carried out for the second- 
order (i.e. two-point) correlation function of the scalar field: there are two 
reasons for this choice. First of all, the second-order correlation function is 
the Fourier transform of the spectrum of scalar variance, a statistical indicator 
widely used to characterize most of the statistical properties of scalar turbu- 
lence. Moreover, for the Kraichnan model only the second-order correlation 
function has a simple, closed analytical expression; for higher-order correlations 
only perturbative expressions (for example in the limit of small £) are available 
|23 El- This latter fact suggests that, while a constant-eddy-diffusivity may 
be found also for them, the cxportability of our refined closure seems ruled out, 
in view of the considerations reported before equation l|37|l . Furthermore, even 
if the functional form of the closure was preserved, a modification of the effec- 
tive coefficients would mean that strong small-scale fluctuations — associated 
to higher-order correlation functions — must be described by parameters dif- 
ferent from the ones used for less intense fluctuations. That would question the 
applicability of closure models to the description of the statistics of turbulent 
fields as temperature or concentration, which are characterized by a wide range 
of fluctuation intensities. 

Let us also mention another possible extension, in the presence of inhomo- 
geneities or anisotropics. Starting from some known results about the two-point 
equal-time correlation function in these nonideal situations E| > it would be 
of interest to derive the corresponding LES closures from first principles. 
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The two above-mentioned closures, and are 'universal' in the sense 
that they were derived from first principles. However, in the remaining subsec- 
tions of we showed that some other closures may emerge. A slight diversi- 
fication in the filtering length between velocity and scalar allows to capture the 
0(l/r) 2 correction even imposing a constant-eddy-diffusivity parameterization. 
On the other hand, it is possible to obtain the same accuracy doing without any 
renormalized diffusion coefficient. Finally, we showed that our refined model 
can be applied to other quantities than the coarse-grained correlation function: 
e.g., with a simple change in the numerical coefficients, one can reproduce the 
original, fully-resolved correlation, i.e. the 'pure' inertial-range scaling, even in 
the presence of finite-size effects. 

A last question that naturally arises is whether these results are relevant for 
realistic advection models, e.g. relaxing the assumption of ^-correlation in time. 
Of course, the value of the effective diffusivity and of the triplet of numerical 
coefficients that appear in these closures can be analytically computed only in 
the Kraichnan model. However, we believe that the form of the parameterization 
can be exported without modifications to real situations as well. Clearly, in this 
case, such free parameters have to be determined a posteriori by some empirical 
procedure. The validity of this approach has already been checked by direct 
numerical simulations |28|. which confirm such a preview even with velocity 
fields that are real solutions of the NS equation. 
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A Exact analytical expressions 

In this section we show the exact analytical expressions of some quantities, of 
which we have only reported perturbative expressions throughout the paper. 
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A.l Coarse-grained passive scalar 

The correlation of the coarse-grained passive scalar mentioned in (|21|) is: 
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where R± = r ± 21 and i? = sgn(i?_). The fuse-point value is given by 

48F l 2 -t 



2«(2- 0(4 -0(5 -0(6 -0^o ' 



i.e. (6 2 ) < (9 2 ) = c. 

Keeping lfl%|l and f2T|l into account, this leads to the following expression for 
the two-point coarse-grained scalar structure function in the inertial range: 
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(40) 



Expression Ij40|l shows that spurious corrections 0(l/r) 2 ~^, which are absent 
in the correlation, appear in the structure function when filtering. In other 
words, it is not true that the coarse-grained structure function is obtained by 
performing a double convolution on the fully-resolved corresponding quantity 
(this result is in accordance with the integral expression of the coarse-grained 
velocity structure function in \22l ). Of course, this fact does not spoil our 
closures at all, but it makes it easier to deal with correlation functions. 



A. 2 Coarse-grained velocity 

The coefficients A(r) and B(r) appearing in (|23|l are uniquely determined by 
the system 

f 2rB'(r) - 6B(r) = -rG{r) 
\ 3A(r) + B{r) = G{r) 
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with 5(0) = and 
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A series expansion in l/r gives 



A(r) = D Q r^ 
and 



(2 + 0- 



B(r) = D r C 
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B Limit cases: purely-diffusive or smooth flow 

The values £ = and £ = 2 have been excluded from our analysis up to now, 
because they represent two limit cases: a purely diffusive flow and a smooth 
flow, respectively. It is however interesting to analyse them, as in the former 
case an exact closure can be found, and in the latter a logarithmic law arises. 

B.l The case f = 

In the case £ = the definition of 77 (|15[) is meaningless, but physically it 
corresponds to a diffusive range extending to infinity. This can be understood 
very simply noticing that the velocity field reduces to a white noise, as its 
second-order moment (I14f) takes a diagonal form completely independent from 
the separation: Dj?J(r) = (d — 1)D Q S^. 

Consequently, the overall effect of the advective term is barely of diffusive type 
and consists in an addition of the quantity {d— 1)Dq to the original contribution 
2kq in (|16|l . As there is no more need to split the integration interval in (|17(l . 
the correlation is exactly expressed by 



Cf (r) 
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for r < L 
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for d ^ 2 



19 



and 



F 



2d[2n + (d- 1)D ] 
By definition, an exact calculation in 3-D yields: 

2" 

C?'(r) =c-kr 2 ' 



with 



kr 2 



10(/s + A>) ' ' 

In this case, the structure functions of the coarse-grained field and of the origi- 
nal one turn out to be identical: S% \r) = 2kr 2 = s!f\r). 
The interesting point is that CED closure (|29ll is now exact also for the cor- 
relation function, as it happened in (|llf) only for the mean value. This is in 
accordance with the vanishing of the last term on the right-hand side of (|36|) . 
From the analytical point of view, the exactness of the eddy-diffusivity closure 
is due to the vanishing of the second-order structure function of the coarse- 
grained velocity, D^J (r) = 0, which imposes an exact balancing between the 
forcing and the diffusive terms in (|32ll . 



B.2 The case f = 2 

If £ = 2, the second-order spatial increments of the velocity (|14|) scale with r 2 , 
so v is a differentiable field. The passive-scalar correlation is exactly given by 

d£ \r) = c- fcln T —„ + \ for r < L , (41) 



where 



and 



L 2 + T] 2 

Fo 
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d 2 {d-l)Do 2 l \^2^ ' 2 : (d-l)D L 2 / 
2F1 being the hypergeometric function. 

Exploiting its definition, for d — 3, the coarse-grained scalar correlation shows 
a second-order correction in l/r with respect to its fully-resolved corresponding 
value i|4*TJl . 

cf(r)=C<«- JrfzV + ofrV- (42) 



d d 2Ka 



30D 

while no difference exists for the velocity structure function: D^J — D^) . 
The key point is that 0(l/r)£ and 0(1 /r) 2 obviously coincide, so CED closure 
cannot be introduced by itself because it is intimately entangled with the im- 
proved closure (|3"S|l . The latter captures the second-order correction in ll4*2l 
correctly. 
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